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Non-perturbative solutions to the quantum-field theory is a topic of current and broad interest, 
especially for the heavy ion and laser physics communities, since they investigate particle pro- 
duction in the presence of strong external fields. We have solved a non-perturbative QED + 
external field problem of an electron in a strong transverse confining potential using Hamiltonian 
light-front quantum field theory in a basis function representation. The invariant mass spectra and 
the anomalous magnetic moment of the lowest state for this two-scale system are also evaluated. 
With this method the perturbative QED results are also reproduced with a good accuracy. We also 
discuss the extension of the method to other problems as well. 
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1. Motivation 



The Hamiltonian light-front formalism is ideal for solving problems non-perturbatively, since 
time is set along the light-front and evaluated experimental observbles, such as masses, form fac- 
tors, and structure function, are Lorentz frame independent. Light-front Hamiltonian quantum field 
theory also has similarities with non-relativistic quantum many-body theory. This connection was 
exploited in [JIJ], where a "Basis Light Front Quantized (BLFQ)" approach was outlined by adopting 
a light-front single-particle basis space consisting of the 2-D harmonic oscillator for the transverse 
modes and a discretized momentum space basis for the longitudinal modes. Adoption of this basis 
is also consistent with recent developments in AdS/CFT correspondence with QCD 

In [||] this approach was applied to address the problem of an electron in a transverse harmonic 
cavity with the QED Hamiltonian evaluated on the light-front in a Fock space consisting of electron 
states and electron plus photon states. The external field was included non-perturbatively and the 
eigenvalues, eigenvectors and anomalous magnetic moment were solved. In a non-renormalized 
case, the obtained electron anomalous magnetic moment was within 1.5% of the theoretically ex- 
pected result (Schwinger moment), when extrapolated to the zero external field limit. Applying 
a sector-dependent renormalization scheme [||] to this Hamiltonian, the zero external field results 
were consistent with related works Refs.|^, |6j] and Refs.ffT], §], where the one-photon truncated 
light-front Hamiltonian was regulated with a Pauli-Villars regularization scheme. 

The nonperturbative analysis presented in [^] could be applicable to measurements of the 
(gyromagnetic) ratio of the spin precession to Larmor frequencies of a trapped electron in strong 
external electromagnetic fields, and can be straightforwardly extended by incorporating higher 
Fock-space sectors. It also serves as a first step towards the studies of non-perturbative QED 
relevant for the anomalous enhancement of lepton production at RHIC [[J and for proposals for 
producing super-critical fields with next-generation laser facilities fliO| , O]. Another direction for 
QED applications follows the lines of [12, 13], where Fourier transform of the Deeply Virtual 
Compton Scattering amplitude with respect to the skewness variable £ at fixed invariant momentum 
transfer was observed to be analogous to the diffractive scattering of a wave in optics. In analogy 
with this "hadron optics", the light-front electron wave functions computed in [Q] can be used to 
evaluate the form factors of the electron and thereby introduce "electron optics". 

Most importantly, the research in ||] also serves as a foundation for solving other quantum 
field theories at strong coupling, such as the light-front QCD Hamiltonian in the nonperturbative 
domain. In order to extend the research to QCD, methods for treating the color degree of freedom 
in a computationally efficient manner were already introduced and evaluated in [[!]] . The next step 
is then to incorporate the color degree of freedom into the quantum field theory code, that has 
already passed an important accuracy test in 

We will next review the method used in [[T], |[], adding some details not previously discussed. 



2. BLFQ Hamiltonian Framework 

In BLFQ approach the Hamiltonian is expressed in terms of basis functions, and the size of the 
resulting Hamiltonian matrix is regulated by imposing both physical symmetries and different cut- 
off conditions for the basis states. Increasing the size of the basis will inevitably lead to substantial 
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computational requirements both in the computation of the matrix elements themselves and in the 
diagonalization of the matrix. In order to be able to extrapolate the results to the continuum limit, 
where the cut-offs are removed, a rapid convergence of the results is then highly desirable. 

We define our light-front coordinates as x ± = x° ±x 3 , x^ = (x l ,x 2 ), where the variable x + 
is light-front time and x~ is the longitudinal coordinate. We adopt x + = 0, the "null plane", for 
our quantization surface. In our choice of framework [p], we quantize QED on the light-front 
using the light-front gauge, and add the harmonic oscillator potential in the transverse direction to 
confine the system in those directions. To simplify the numerical work, we choose the transverse 
basis function scale and the trap scale to coincide. As a consequence, we cannot obtain zero 
external field QED results directly, but via extrapolation, as shown later. 

Our basis states consist of 2-D harmonic oscillator (HO) states, which are combined with 
discretized longitudinal modes, plane waves satisfying selected boundary conditions. The HO 
states are characterized by a principal quantum number n, orbital quantum number m and HO 
energy O. We express the 2-D oscillator in momentum space as a function of the dimensionless 
variable p = \p L \ / \J MqQ., where Mq has units of mass. The orthonormalized HO wavefunctions in 
polar coordinates (p, <p) are then given in terms of the Generalized Laguerre Polynomials, ij, (p 2 ), 
by 



^nmip^ =*nm(p,<P) = (P?M = yj ^ J V" <' "' J C (p 2 ) , (2-D 

with eigenvalues E n m = (2n + \m\ + 1)£2. These wavefunctions are orthogonal and form a complete 
set of states, thus 

£<J>: m (^)0„ m (^) = (2tt) 2 S( 2 V (2-2) 

nm 

The longitudinal modes, Yk, in our basis are defined for — L <x<L with periodic boundary 
conditions (PBC) for the photon and antiperiodic boundary conditions (APBC) for the electron: 

Vk (x-) = ^=^ kx \ (2-3) 
V 2L 



where k = 1,2,3, ... for PBC (we neglect the zero mode) and k = |, §, §, ... for APBC. The full 3-D 
single-particle basis state is defined by the product form 

^,„, m (^,P,<P) = Vk{x~)®nAP,<?)- (2-4) 



Following Ref.[14] we introduce the total invariant mass-squared M for the low-lying phys- 
ical states in terms of a Hamiltonian H times a dimensionless integer for the total light-front mo- 
mentum K 

M 2 + PjPjl -> M 2 + const = P + P = KH (2.5) 

where we absorb the constant into M 2 . The non-interacting Hamiltonian Hq for this system (where 
now the transverse functions for both the fermion and the boson are taken as eigenmodes of the 
trap) is then defined by the sum of the occupied modes i in each many-parton state as 

hq = 2m p- = 2 -^y + N + 1 + (26) 

K ; Xj 
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where m,- is the mass of the partem i. The photon mass is always set to zero in the following and the 
electron mass m e is set at the physical mass 0.51 1 MeV in our non-renormalized calculations. We 
also set Mo = m e . 

The basis is limited to fermion and fermion-boson states, so the light-front QED Hamiltonian 
interaction terms we need are the fermion to fermion-boson vertex, given as 



V e ^ey = 8 J dx + d 2 x ± ^(x)y^(x)A^(x) 



(2-7) 

X+=0 



and the instantaneous fermion-boson interaction, 

,2 



Vey^ey = y J dx + d 2 X ± ^ (f 'A v *¥) 



(2.8) 

x+=0 



where the coupling constant g 2 = 4na, and a is the fine structure constant taken to be a = l37 \ m - 
When expressing the free fermion and boson fields in terms of our basis functions, the complete 
set of quantum numbers needed to specify a state are a = (k,n,m,X), where A is the helicity. The 
fermion and boson annihilation operators are then written as 

b(p^,a) = Y,b(oc)^nm(p^), (2.9) 

nm 

a(p X ,Y)=Z a (y^nm{p X ), (2-10) 

nm 

where 

[a(f),aHy>)] = tf, (2.1 D 

{b(a),b\oc')} = 8i', (2.12) 

and the truncated set of quantum numbers a = (k,X). After this replacement and proper nor- 
malization, the non-spinflip vertex terms of Eq.( |2.7[ ) are oc MqQ., similar to the non-interacting 
Hamiltonian of Eq.(^6|), whereas spinflip terms are oc ^MoCbn e . Selecting the initial state fermion 
helicity in the single fermion sector always as "up", the process e — > ey is nonzero for 3 out of 8 
helicity combinations, and the process ey — > <?yis nonzero only when all 4 spin projections aligned 
(2 out of 16 combinations). The resulting Hamiltonian matrix is thus sparse. 

As a symmetry constraint for the basis we fix the total angular momentum projection J z = 
M + S = where M = Y*i m i i s th e total azimuthal quantum number, and S = tne total spin 
projection along the x~ direction. For cutoffs, we select the total light-front momentum, K, and the 
maximum total quanta allowed in the transverse mode of each one or two-parton state, N max , such 
that 

£2n ; -+K-| + l<AU x , (2.14) 



where, for example, kj defines the longitudinal modes of Eq.(|2.3[) for the f" parton. Eq.(|2.13|) 



signifies total light-front momentum conservation written in terms of boost-invariant momentum 
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fractions, x,. Since each parton carries at least one unit of longitudinal momentum, the basis is 
limited to K partons. Furthermore, since each parton carries at least one oscillator quanta for 
transverse motion, the basis is also limited to N max partons. Thus the combined limit on the number 
of partons would be mm.(K,N max ), if the Fock space was not truncated. 



3. BLFQ Hamiltonian results 

Here we present some numerical results from [^J for cases where the cutoffs for the basis 
space dimensions are selected such that K increases simultaneously with the N max . The resulting 
dimension of the Hamiltonian matrix increases rapidly. For N max = K = 2, 10,20, the dimensions 
of the corresponding symmetric d x d matrices are d = 2, 1670, 26990, respectively. Since we 
employ a mix of boundary conditions and all states have half-integer total K, we quote /f-values 
rounded downwards for convenience, except when the precise value is required. 
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Figure 1: Eigenvalues (multiplied by K) for a non-renormalized light-front QED Hamiltonian which in- 
cludes the fermion-boson vertex and the instantaneous fermion-boson interaction without counterterms (fig- 
ure adapted from Qj). Closeup of the lowest-lying eigenvalues on the right panel. The basis is limited to 
fermion and fermion-boson states satisfying the symmetries. The cutoffs for the basis space dimensions are 
selected such that K increases simultaneously with the Nmax- 

In the left panel of Figjl] we show the eigenvalues (multiplied by K) for a non-renormalized 



light-front QED Hamiltonian given in Eqs.( |2.6| , |2.7| , |2.8| ), with fixed Q. = 0.05 MeV. In the right 



panel we show a closeup of the lowest-lying eigenvalues only. These eigenvalues correspond to 
a solution dominated by the electron with n = m = 0, and are expected to be KEq = ml+MoQ., 
where the latter term accounts for the lowest state of transverse motion of the electron allowed in 
the chosen basis. In Fig.[T] the lowest eigenvalues for a fixed £2 fall below that value, as the size 
of the matrix increases. The contribution of the spinflip terms to the lowest eigenvalues is very 
small, and, as a result, the lowest eigenvalues for a fixed N max = K depend linearly on Q.. Since our 
system is in an external field, the lowest physical mass eigenstate (not known experimentally) can 
deviate from the free-space mass. Therefore, before renormalization, we choose only to consider 
cases where the mass eigenvalue falls within 25% of the free electron mass. 
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The ordering of excited states in Fig.[I], due to significant interaction mixing, does not always 



follow the highly degenerate unperturbed spectrum of Eq.p.q). States dominated by spin-flipped 
electron-photon components are evident in the solutions. Nevertheless, the lowest-lying eigen- 
values appear with nearly harmonic separations in Fig.[j] as would be expected at the coupling of 
QED. The multiplicity of the higher eigenstates increases rapidly with increasing N max = K and 
the states exhibit stronger mixing with other states than the lowest-lying states. In principle the 
fermion-boson basis states interact directly with each other in leading order through the instanta- 
neous fermion-boson interaction, but numerically the effect of this interaction is very weak, and 
thus does not contribute significantly to the mixing. Even though we work within a Fock space 
approach, our numerical results should approximately equal the lowest order perturbative QED 
results for sufficiently weak external field. 
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Figure 2: (Color online) Square root of the (scaled) election anomalous magnetic moment as a function of 
the ttansverse external field for a sequence of increasing basis spaces indicated in the legend. These are non- 
renormalized results where the mass eigenvalue falls within 25% of the free electron mass. Extrapolation 
to zero external field yields 0.1 121, compared with the theoretical 1-loop QED prediction ("Schwinger") of 
0.1 125. Figure is adapted from [^J. 

In Fig.^ we show the results for the square root of the electron anomalous magnetic moment 
(scaled), \J 8pL /g 2 , as a function of Q. obtained from the lowest mass eigenstate. That is, we plot 
the magnitude of the probability amplitude that the electron has its spin flipped relative to the single 
electron Fock-space component in the range where the results are converged. For even N max = K 
the results converge rapidly for N max = K > 14. The results for odd cutoffs (not shown) track even 
cutoff results as N max = K increases. Below Q. < 0.05 MeV (results not shown), in the weak external 
field region, all the interactions are quenched at fixed N max = K, and not converged, due in part to 
our requirement that the HO basis states track the external field. In order to compare our results 
with the square root of the ratio of the Schwinger result ^ for the coupling constant g 2 = 4na, we 
perform an extrapolation of the above results for N max = K = 12, ... ,20 to the zero external field 
limit £2 = Mev. An excellent agreement with the results is obtained by a fit function = 
a(l + M2 2 + c£2 4 )exp(— d£i), with a = 0.1121. This is < 1% deviation from the Schwinger result 
of 0.1 125, which is reasonable in light of our numerical accuracy and extrapolation uncertainties. 
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Figure 3: (Color online) Individual fits to the renormalized results for square root of the (scaled) electron 
anomalous magnetic moment for N max = K = 10, ... , 20. The inset shows the continuum limit extrapolation 
of the zero external field results from the main panel as a function of 1 /K. Figure is adapted from 



In Fig .[3] we renormalize our results for a/ dj^/g 2 by applying a sector-dependent normalization 
scheme from Ref.[Q|. In our present limited Fock space, we need only the mass counterterm 8m e . 
This 8m e is added to the mass term in the diagonal one-electron part of the Hamiltonian Eq.(^). In 
the absence of a known experimental mass for renormalization due to the external field, we adjust 
8m e such that the lowest eigenstate remains at KEq = m 2 . + MqCI. To eliminate possible effects 
from the quenched interactions at small external fields, we only include results with the external 
field Q. > 0.2 MeV. Again, individual fits of the form given Fig.|| in are an excellent representation 
of our results. The range of the extrapolated values is 0.1077 < a < 0.1216. 

The convergence with an increasing cutoff is now less rapid than in the non-renormalized case 
shown Fig.|[ In order to approach the continuum limit N max = K — > °°, we perform further extrap- 
olation to the zero-£2 results of Fig.|3| The inset of Fig.|| shows linear extrapolation of the results 
of the main figure in l/K to the continuum limit N max = K — > °°. To verify the stability of the 
results, an extrapolation based on the £2 > 0.1 MeV fits (not shown) is also given. The extrapolated 
continuum values are 0.1362 (0.1383) for £2 > 0.2(0.1), respectively, and thus about 20% above 
the Schwinger result 0.1125. As mentioned in Sec.|], an enhancement of this magnitude was also 
observed in related works, Ref.[|]] and Refs.^J, |8|], where one-photon truncated light-front Hamil- 
tonian was regulated with Pauli-Villars (PV) regularization scheme. With PV regularization as well 
as in our renormalized results, interpreted from a perturbation theory perspective, the intermedi- 
ate state propogators are developed from a dynamical (non-perturbative) electron mass rather than 
using the unperturbed mass needed for direct comparison with perturbation theory. Thus one may 
appreciate why the renormalized results are distinct from the lowest order perturbative results. 

The extension of this method both to a larger Fock-space and to QCD will proceed as outlined 
in [[I]]. In the case of non-interacting QED, the state density as a function of the state energy E, given 
by Eq.( |2.6| ), was found to increase exponentially as N max = K increases. Similarly, implementation 
of the color dramatically increases the state density over the case of QED, as shown in Fig.|], but it 



7 



Electron in a transverse harmonic potential 



H. Honkanen 



was also found, that use of a global color-singlet constraint is effective in minimizing the explosion 
in basis space states. We anticipate to be able to handle both of these challenges with the parallel 
codes developed, tested and applied in [p], ||], that compute the Hamiltonian matrix and solve for 
its eigenvalues and eigenvectors 




Figure 4: Number of color space states that apply to space-spin configuration of selected multi-parton 
states for two methods of enumerating the color basis states. The upper curve shows counts of all color 
configurations with zero color projection. The lower curve counts global color singlets. Figure adapted 
Ref.[|. 



References 

[1] J. P. Vary et ai, Phys. Rev. C 81,(2010) 035205. 

[2] G. F. de Teramond and S. J. Brodsky, Phys. Rev. Lett. 102, 081601 (2009). 

[3] H. Honkanen, P. Maris, J. P. Vary and S. J. Brodsky, arXiv: 1008.0068 [hep-ph]. 

[4] V. A. Karmanov, J. F. Mathiot and A. V. Smirnov, Phys. Rev. D 77 (2008) 085028. 

[5] S. J. Brodsky, J. R. Hiller and G. McCartor, Phys. Rev. D 58, 025005 (1998). 

[6] S. J. Brodsky, V. A. Franke, J. R. Hiller, G. McCartor, S. A. Paston and E. V. Prokhvatilov, Nucl. 
Phys. B 703 (2004) 333. 

[7] S. S. Chabysheva and J. R. Hiller, arXiv:091 1.3686 [hep-ph]. 

[8] S. S. Chabysheva and J. R. Hiller, Phys. Rev. D 81 (2010) 074030. 

[9] A. Adare et al. [PHENIX Collaboration], Phys. Rev. C 81 (2010) 03491 1. 

[10] M. Ruf, G. R. Mocken, C. Muller, K. Z. Hatsagortsyan and C. H. Keitel, Phys. Rev. Lett. 102 (2009) 
080402. 

[11] C. K. Dumlu and G. V. Dunne, Phys. Rev. Lett. 104 (2010) 250402, and references therein. 

[12] S. J. Brodsky, D. Chakrabarti, A. Harindranath, A. Mukherjee and J. P. Vary, Phys. Lett. B 641 (2006) 
440 [arXiv:hep-ph/0604262]. 

[13] S. J. Brodsky, D. Chakrabarti, A. Harindranath, A. Mukherjee and J. P. Vary, Phys. Rev. D 75 (2007) 
014003 [arXiv:hep-ph/0611159]. 

[14] S. J. Brodsky, H. C. Pauli and S. S. Pinsky, Phys. Rept. 301 (1998) 299. 



8 



